Random Delaimay triangulations 
and metric uniformization 

Gregory Leibon 
February 8, 2008 



1 Introduction 

The primary goal of this paper is to develop a new connection between the 
discrete conformal geometry problem of disk pattern construction and the 
continuous conformal geometry problem of metric uniformization. In a nut- 
shell, we discuss how to construct disk patterns by optimizing an objective 
function, which turns out to be intimately related to hyperbolic volume. 
With the use of random Delaunay triangulations we then average this objec- 
tive function to construct an objective function on the metrics conformal to 
a fixed one. Finally using this averaged objective function we may reprove 
the uniformization theorem in two dimensions. 

This paper is organized as follows. In section |2| we introduce the use of 
random Delaunay triangulations by presenting a new proof of the Gauss- 
Bonnet formula. In section |3| we present the disk pattern ideas and the 
objective function in the disk pattern setting. In section [| we use the tech- 
niques developed in section |2| to average the objective function from section 
[| and produce the objective function on metrics. In section || we point to 
some open questions. 

This article is based on the author's thesis ||, where readers can find a 
detailed treatment of everything that takes place here. 
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2 A Random Proof of the Gauss-Bonnet The- 



orem 



Recall that the Gauss-Bonnet formula may be stated for a compact bound- 
ary less Riemannian surface M as: 



with k the Gaussian curvature and x(-M) the surface's Euler characteristic. 
The Euler characteristic is a topological invariant which can be computed 
relative to any triangulation via 



with F, E, and V the number of faces, edges and vertices in the triangulation. 
The proof here is accomplished by randomly triangulating the surface and 
then noting that while x(M) is constant F, E, and V are now random vari- 
ables and have expected values which can be computed. As the density of the 
randomly distributed vertices goes to infinity one finds these expected values 
produce the Gauss-Bonnet formula, along with a probabilistic interpretation 
of Gaussian curvature. 

To begin with, we need to define what we mean by a random triangulation 
of a surface. The first step is to ignore the fact that anything random is going 
on here and to simply attempt to construct a geodesic triangulation in a fixed 
metric g from a given set of points, p = {pi, ...,p n }. To accomplish this we 
produce an abstract two complex by examining all the triples and pairs in 
{pi, ...,p n } and deciding whether or not to put in a face for a given triple or 
an edge for a given pair. This decision procedure will be relative to a certain 
positive number 5 - the decision radius. The procedure is to put in a face 
for a triple or an edge for a pair if the triple or pair lies on a disk of radius 
< 5 which has its interior empty of points in {pi, ...,p n }. If a triple is on an 
empty disk then each pair in this triple is on an empty disk, so we indeed 
have a 2-complex. This proceedure is called Delaunay's "empty sphere" 
method and was introdueced in ||. It is elementary to see that there is a 
positive decision radius such that one can view the edges of this 2-complex 
as geodesies in the g metric. When this procedure forms a triangulation of 
M we call the resulting triangulation the Delaunay triangulation. 

At this point it is useful to introduce a geometric criterion on a set of 
points {pi, ...,p n } guaranteeing that it produces a Delaunay triangulation. 




X {M) = F-E + V 
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Definition 2.1 We will call a set of points {pi, ...,p n } generically 5-dense if 
each open ball of radius 5 contains at least one Pi and {p%, ...,p n } contains 
no four points on a disk of radius less than 5. 

It is straightforward to see that 

Lemma 2.2 There is a 5 > such that if {pi, ...,p n } is generically J-dense, 
then {pi, ...,p n } forms a Delaunay triangulation. 

Now enters the randomness. From lemma |2.2| , when our points are dis- 
tributed with a high density we expect they will typically form Delaunay 
triangulations. This procedure is somewhat independent of our choice of 
distribution, but for all that takes place here we will assume that we are 
using a Poisson point process relative to a density denoted A (see remark 
j2.8| ). Letting E A (L) denote the expected value of a random variable L and 
letting 0(X~°°) mean a quantity decaying faster than any polynomial in A, 
we indeed have that... 

Lemma 2.3 The probability that a set of points form a triangulation is 
1 + 0(A~°°). Also, if L is any one of the random variables E, V, or F, then 
L has expected value equal to E A + 0(A~°°), with L the random variable 
which is L when the points form a triangulation and zero otherwise. 

From this lemma one sees that the constant x(M) satisfies 

X (M) = E A (xW)) + ©(A" 00 ) = E A (V - E + F) + 0{\-°°). 
Applying the fact that expected values add gives us 

X (M) = E A (F) - E A (E) + E A (T) + 0(A-°°). 
Furthermore, in an actual triangulation we have \E = F so 

X (M) = E A (F) - ^E A (T) + O(A-). 

Since A is the density, if we let A denote the area of M then the expected 
number of vertices is XA. This observation along with lemma |2.3| gives us 

X (M) = A\-±E x (F) + 0(\-°°). 

So we have... 
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Formula 2.4 (Euler-Delaunay-Poisson Formula) 

X (M) = lim (a\ - \-Ex{F) 

The goal now becomes to compute E\(F). Let 5 be small enough to 
satisfy lemma |2[^, let Vg C x 3 M = M x M x M he the set of ordered triples 



living on circles of radius less than 5, and let a(y) be the area of the disk 
associated to the triple y G Vs. If you are familiar with how to compute using 
the Poisson process you will find 

E x (F) = l ( e' Xa ^(XdA) 3 . (1) 
6 Jv 6 

(If you are not familiar with this situation you should see remark 2~1L) 

In order to explicitly compute the integral in formula |l| it is necessary 
to put coordinates on Vs. To accomplish this, first one chooses a way to 
discuss directions at all but a finite number of tangent planes of M (via an 
orthonormal frame field). Then one can parameterize a full measure subset 
of V s with a subset of {(#i, 2 , 3 , r,p) E S 1 x S 1 x S 1 x (0,5) x M} by 
starting at the point p G M and moving a distance r in each of the three 
directions 9\, 6 2 , and 63 . Note that when fixing p and r and varying we 
produce a Jacobi field, whose norm we shall denote jV Using this notation, 
letting d9 = dOidd^ddz, and letting v{9) be the area of the triangle in the 
Euclidean unit circle with vertices at the points corresponding to the {Oi}, a 
straightforward computation shows us that formula [I] is 

%f f I e- x ^jeMeXO)dedrdA. (2) 

O J M JO Jx^S 1 

in these coordinates. 

The Taylor expansions of a(9, r,p) and jg are controlled by the Gaussian 
curvature up to the fourth and third order terms respectively. So, after 
potentially shrinking 5 a bit to exploit this control, we may Taylor expand, 
integrate, and apply the mean value theorem to express as 

E A (F) = 2AX — — [ kdA + o(\-$). (3) 

IX J m v ' 



In particular, formula |2.4j may now be plugged into the Euler-Delaunay- 
Poisson Formula to give simultaneously a probabilistic interpretation of cur- 
vature and a proof of the Gauss-Bonnet theorem. 
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Theorem 2.5 (Euler-Gauss-Bonnet-Delaunay Formula) 

X (M) = lim (AX - ^E A (F)) = i / fcdA. 

Note this allows us to interpret the Gaussian curvature as the density of 
the defect in the expected number of faces in a random Delaunay triangu- 
lation in the surface's geometry relative to what would be expected in the 
Euclidean plane. 

Remark 2.6 It should be noted that the lemmas and computations above 
are completely elementary with the exception of the use of the following bit 
of geometry. 

Lemma 2.7 (The Small Circle Intersection Lemma) There is a 5 > 
such that if a triple of points lies on the boundary of a disk with radius less 
than 5, then this disk is unique among disks of radius less than S. 

This lemma is at the core of the proof of lemma |2^, as well the justification 
for the well-defmedness of a(y) and the parameterization in formula ||]. A 
couple of proofs of this lemma can be found in ||, where it is shown that 
we can get an explicit grip on the necessary S by using any 5 < mm{|,r}, 
where i is the surface's injectivity radius and r the surface's strong convexity 
radius. 

Remark 2.8 Any reasonable form of point distribution will fare as well as 
the Poisson point process with regards to this entire proof. This includes 
the uniform distribution of n points, which is simply the product measure on 
M x . . .x M — x n M. This uniform choice in fact eliminates all but the most 
basic probabilistic thinking, though not without a certain aesthetic sacrifice. 

Remark 2.9 Here we will give an informal description of the Poisson point 
process and the derivation of formula [1]. To describe the Poisson point pro- 
cess, imagine breaking the surface up into pieces of size dA small enough so 
that the probability that one of these pieces contains more than one point 
is negligible. Denote one of these little chunks by q and let X q be 1 if the 
region q contains a point and otherwise. The heart of the Poisson point 
process is the assumption that the X q are independent and the probability 
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that q contains a point is E\(X q ) = XdA. Exploiting this independence we 
see that the probability of a region R of area A being empty of points is 



Let R t be the function which is one if the "disk" formed by a "triple" t = 
ill, 12, Q3} is empty of other points, and zero otherwise (the quotations are 
used because it is only after taking the limit that we really have a triple of 
points and an actual disk). In particular note that R t has expected value 
e --W*) Hence by independence we have that 



is the probability of a "triple" forming a face. In particular by the linearity 
of expected values 



which gives us formula [3] in the limit. 

3 Discrete Conformal Uniformization 

The goal here is to describe the needed disk pattern ideas. In order to 
facillitate the use of these ideas in the metric world we will present our 
discrete objects as natural data living within actual geodesic triangulations of 
Riemannian surfaces. For example, imagine starting with a triangulation of 
a hyperbolic surface, by which we mean a geometric surface having constant 
Gaussian curvature — 1. Let s and t be two triangles in this triangulation 
which share an edge e, and denote the complement of the intersection angle 
between the disks in which s and t are inscribed as ip e . See Figure 1. The 
key is to note that ip e can be written down in terms of the angles within the 
triangles as 

Formula 3.1 



E x (U q£R X g ) « (1 - XdA) 



A \A 
dA Ki e . 



E\(R{ qu q2 )q3 }X qi X q2 X q3 ) = e Aa(t) (AcL4) 3 



E A (F)= £ E A (J2 t )= £ e a{t \\dA)\ 




with 



2 



using the notation in the first figure. 
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Figure 1: The notation. 



Comment on Proof: This formula always holds on a constant curvature 
surface, and a proof can be found in ||. For our purposes it is useful to 
introduce the key object showing up in the proof of the negative curvature 
case. This object is an ideal hyperbolic prism. The prism is constructed 
from the angle data of a hyperbolic triangle, namely a set of positive angles 
{A, B, C} such that A+B+C— n < 0. To construct it first form a hyperbolic 
triangle with the {^4, B, C} data, then place the triangle on a hyperbolic plane 
and then place the plane in hyperbolic three space. Now union this triangle 
with the geodesies perpendicular to this 2-plane going through the vertices 
of the triangle. The prism of interest is the convex hull of this arrangement. 
See Figure 2. 

Now suppose we have any triangulation of a Riemannian surface with 
varying curvature; What data may we collect? Well, we may collect the 
topological triangulation T and all the triangle angles. From formula |3T 
this triangle angle data may be organized into the collection of partial angles 

whose values we will identify with the coordinates in a 3F dimensional 
real vector. It should be pointed out that the angles inside the triangles are 



linearly determined by the tpt v ^ a inverting formula 3.1. To be the data of 



a Delaunay triangulation forces certain natural linear constraints, which we 
will capture with the following definition. 
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Figure 2: The prism. 



Definition 3.2 Let an angle sytem be any point x in the above vector space 
such that rr's associated triangle angles are in (0, 7r) and sum up to 2n at any 
vertex. A Delaunay angle system is an angle system x where the complement 
of the informal intersection angle, vf) e (x) = ip*(x)+ip%(x), has values in (0, n). 

We call two angle systems x and y conformally equivalent if they share 
the same informal intersection angles. We record this with idea the following 
definition. 

Definition 3.3 x and y will be called conformally equivalent if ip e (x) = 
i/j e (y) for all e. 

Motivated by the Gauss-Bonnet formula we also define the following. 

Definition 3.4 The curvature of a triangle t relative to an angle system x is 
defined to be 7r subtracted from sum of the triangle angles in t determined by 
x. If every triangle has negative curvature relative to x, then x will be said 
to have negative curvature. The set of Delaunay angle systems with negative 
curvature will be denoted N x . 

There is a simple set of linear equations which will gaurentee that a 
Delaunay angle system is conformally equivalent an angle system in N x . Let 
S be a set of triangles in T, and denote the cardinality of S as \S\. 
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Definition 3.5 An angle system x will be called teleportable if for any set 
of triangles S we have 

5>-^(z))>7r|s|. 

eeS 

As an immediate consequence of theorem 2 in |J we have the following 
lemma. 

Lemma 3.6 (The Discrete Teleportation Lemma) A Delaunay angle 
system x is teleportable if and only if x is conformally equivalent to a negative 
curvature Delaunay angle system. 

The proof of this lemma is completely linear in nature. 

In the setting here we will only be concerned with Delaunay angle sys- 
tems satisfying the conditions of this lemma, which includes the angle data 
associated to a Delaunay triangulation of surface with varying negative cur- 
vature. 

The pleasure derived from the Delaunay angle systems comes from a won- 
drous objective function which lives on "N x . Let Vt(x) denote the volume of 
the ideal hyperbolic prism constructed from a triangle t's angle data rela- 
tive to x, as described in the previous section. Now simply let the objective 
function be 

H(x) = Y,V t (x). 
teT 

The wonder of this function can best be felt by examining its differential. To 
compute this use x's angle data in t to construct a hyperbolic triangle and 
let lf(x) denote the length of the edge e in this triangle. In |§ the following 
formula is produced: 



Formula 3.7 



. , /cosh(/f(x)) - 1\ 
dH = - > log dtp, 




e 



The tangent space at any point of N x is precisely the set of directions 
preserving the condition that the ifj e are constant, hence is spanned by vectors 



in the form C e = -J=^ — So from formula |3~77| we see that a; is a critical 



point of the objective function in a conformal class if and only if 

o = dmo = log t^m^l) _ log fco^O-f 



9 



for all edges e. Thus at a critical point of H we have that lf(x) = l e s {x), 
so the set of hyperbolic triangles formed from x's angle data fit together to 
form an actual hyperbolic surface. 

Definition 3.8 A Delaunay angle system which is the angle data of a hy- 
perbolic surface will be called uniform. 

The question becomes: how many (if any) uniform structures can be asso- 
ciated to a given angle system? The objective function can be analyzed and 
using a compactness argument with boundary control, we find this objective 
function always achieves its maximum (see |§). H is also observed to be 
strictly concave down, so in fact any critical point is i/'s unique maximum, 
and we have 

Theorem 3.9 (The Discrete Uniformization Theorem) If x(M) < 
and a; is a teleportable Delaunay angle system, then x is conformally equiv- 
alent to a unique uniform angle system. 

Remark 3.10 At this point it may be unclear how the above discussion is 
related to a disk pattern problem. Given a geodesic triangulation of a hyper- 
bolic surface, the disk pattern of interest here is the disk pattern produced 
by the circumscribing disks of the triangulation's triangles, which we will 
call the "empty pattern". Notice in the presence of the "empty pattern" 
we may assign to each edge the value of the intersection angle between the 
circumscribing disks of the triangles sharing this edge. A pattern production 
theorem in this setting is an assurance of the existence of an "empty pattern" 
given a topological triangulation and the specification of a sensible intersec- 
tion angle to each edge. To discover what sensible means it is necessary to 
strengthen lemma [3.6| to 

Lemma 3.11 Given any topological triangulation and set of data -0 e G (0, n) 
satisfying both J2eev i> e — 27r and the teleportability condition, there is y G 
N x satisfying ip e (y) = 4> e - 

Hence from theorem |3.9| , under these hypothesis there is an "empty pattern" 
with intersection angles given by 7r — ip e , and we have solved our disk pattern 
problem. 

This pattern problem is equivalent to a generalization of the convex ideal 
case of the Thurston- Andreev theorem when x(M) < (see |J for the details 
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and various generalizations). In the Euclidean case, this extension was car- 
ried out by Bowditch in using techniques similar to Thurston's original 
techniques found in [|TJ|]. The use of an objective function for solving such 
problems was introduced in Colin de Verdiere's (see Question |Q|) , while 



the use of hyperbolic volume as an objective function for producing disk pat- 
terns has its origin in Bragger's beautiful paper H. Hyperbolic volume is 



also used as an objective function in Rivin's [12 



4 Continuous Conformal Uniformization 

The goal now is to bootstrap from the discrete uniformization procedure in 
the previous section to a procedure for producing a conformally equivalent 
uniform metric on a Riemannian surface. Two metrics on M, g and h, are 
conformally equivalent if h = e 2 ^g for a smooth function 0. For metrics by 
uniform structure I will mean a metric with constant curvature. Our goal is 
to prove the classical result... 

Theorem 4.1 (The Metric Uniformization Theorem) Every metric is 
conformally equivalent to a metric of constant curvature, and this metric is 
unique up to scaling. 

We will always be thinking in terms of a fixed background metric called 
g and will label its associated geometric objects like its gradient, Laplacian, 
curvature, norm, area element, or area as V, A, k,\ ■ | , dA or A. For the 
h = e 2 ^g metric we shall denote these objects with an h subscript. 

For starters let us note in the metric world we still have 

Lemma 4.2 (The Metric Teleportation Lemma) Every metric on a sur- 
face is conformally equivalent to a metric with either negative, positive, or 
zero curvature. 

As in the discrete case, this part of the uniformization procedure is completely 
linear and follows at once from the facts that kh = e~ 2< ^(—A(j) + k) and that 
C°°(M) is the L 2 orthogonal direct sum of A(C°°(M)) and the constant 
functions. With this observation in mind, in our x{M) < world we will 
restrict our attention to metrics with strictly negative curvature, and from 
here on out we will assume h has negative curvature everywhere. 
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As described in section |^, relative to a fixed set of vertices we may apply 
Delaunay's "empty sphere" method and produce both the Delaunay triangu- 
lation and its associated "empty disk pattern" (see remark |3.10| ). A confor- 
mal transformation of a metric preserves infinitesimal circles and the angles 
between them, hence for a dense enough set of vertices a conformal transfor- 
mation of a metric nearly preserves the data in the "empty disk pattern" . In 
particular, when we conformally change our metric the angles in the associ- 
ated Delaunay triangulation will change in a manner strongly resembling a 
discrete conformal change, as introduced in the previous section. With this 
observation in mind, if we choose a dense p = {pi, . . . ,p n } and a topological 
triangulation, T, with p as its vertices, then we could measure how close to 
uniform h is with the objective function of the previous section. Specifically, 
we could let 

H h (p) = EK h (p). 
teT 

where V t h (p) is computed using the angle data associated to the triangulation 
viewed in the h metric. This of course means connecting the needed vertices 
of p with h geodesies and measuring the resulting h angles. 

To capture this dense enough set of vertices it is natural to average 
over all sets of vertices with a fixed density and take the limit as the vertex 
density goes to infinity, as we did for the random variable F in section |2|. 
To be explicit, we will distribute points with a density A relative to g's area 
measure, and replace the R t in remark |2.8| with the function that is V t h if the 
"triple's" disk (in g's metric) is empty of points, and zero otherwise. We will 
denote the expected value as E 9 x (H h ), with the superscript g there to remind 
us of the background metric choice. Let 

I°{h)= lim (Hfc — H g ), 

A— >oo 

where the second term, H g , is independent of h and is needed only to nor- 
malize the computation. From this construction we expect that I 9 is an 
objective function capable of uniformizing a negatively curved metric, and 
to confirm this, it is useful to explicitly compute I 9 . 

Theorem 4.3 If g and h are conformally equivalent with h = e 2r ^g then 

P(h) = _ f | V0| 2 + (A<p - k) log(A<£ - k) + k log IJfcldA 
Jm 
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Proof: In performing this computation, we first easily arrive at the analog 
of formula [I], namely 

I 9 (h) = l - I (V h (A l ,A 2 ,A 3 )-V^A h A 2 ,A 3 ))e- x ^(XdAf, (4) 
6 Jv 5 

with V h (A\, A\, A%) the volume of the prism determined by the triangle 
angles {A\, A\, A^} formed using the h metric. At this point theorem |4.3j 
follows from equation [| exactly as equation ^| followed from equation p]: We 
change coordinates, Taylor expand, integrate, and use the mean value the- 
orem, to arrive at our formula for I 9 (h) plus an error term which in this 
case is of the form O (A~5 log(A)J. It should be noted that this procedure, 
which is straightforward when used to produce equation [| becomes consid- 
erably more involved in this setting. First one must Taylor expand A^(r, 9,p) 
in r, which involves solving the boundary value problem determined by the 
geodesic equation up to the second order. With this one may expand V h in 
r, though care is needed since V^'s differential has singularities. The singu- 
larities can be dealt with, and the needed integrals explicitly computed, to 



arrive at theorem [4.3| . The details can be found in |J. q.e.d 

Now we will use I 9 to mimic the proof of the discrete uniformization 
theorem form the previous section here in the metric world. 



Proof of Theorem |4.1| First we confirm that critical points are uniform 
by computing J 9 's differential. Using the notation of theorem [4.3| , it is natural 
to view I 9 as a function on the possible (f> where h = e 2< ^g. With this view 
point we have that the Frechet derivative of I 9 at </> in the direction if) is 

DP(iP) = - f Aipiog\k h \dA. (5) 

Comment. Equation ^] implies that the flow generated by I 9 is the "log 



Ricci" flow, and an alternate proof of theorem |4.1| is to use this use this flow 
as Hamilton used the Ricci flow in M. 

Back to our optimization proof. A straightforward regularity argument 
along with equation |5| for the Frechet derivative assures us that kh is smooth 
at a critical point. From formula |5| at a critical point log \ kh\ is L 2 orthogonal 
to A(C°°(M)). Recalling once again that C°°(M) is the L 2 orthogonal direct 
sum of A(C°°(M)) and the constant functions, we see that log \ k h \ is indeed 
constant. So, in analogy to the discrete metric is critical if and only 

if it is uniform. 
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Using a well chosen space of candidate metrics we find that a compact- 
ness argument with boundary control guarantees the existence of a critical 
point where I 9 achieves its maximum value. For example one could use the 
"xlog + (x)" Orlicz-Sobolov closure (see ||) of 

V = {(p <E C°° \ [ <j>dA = and - A0 + k < 0}, 

JM 

along with some basic functional analysis to arrive at the needed existence 
statement (see |§ for the details). 

Since H g was constructed out of scale invariant angle data, we can only 
hope for uniqueness up to scaling. As in the discrete case, the uniqueness 
follows from the concavity of I 9 . To be more specific, we use the fact that 
all critical points are smooth and that the Frechet Hessian at <fi £ V applied 
to (ip,ip), 

D 2 P(^) = -/ W\ 2 + j^rdA, (6) 
Jm A</> — k 

is strictly negative at a non-zero if) satisfying J M ipdA = 0. 

So we have proved theorem |4.1| by mimicking the discrete case's argu- 



ments, q.e.d 

Notice that we appear to have an infinite number of objective functions, 
one for each metric g. Fortunately any pair of these objective functions differ 
only by a constant. In order to see this it is useful to recall two other functions 
related to metric uniformization, the log(det(A/,,)) and the metric entropy. 
The log of the determinant of the Laplacian had its uniformization properties 



explored by Osgood, Phillips, and Sarnak in [O, while the entropy 



E(h) = - f k h \og\h h \dA h 
Jm 



turned up Hamilton's paper on surface uniformization 0. As a straightfor- 
ward consequence of theorem and Polyakov's formula (see f|TU|) telling us 



that the log(det(Aft,)) with in a conformal class can be expressed as 
1 /l 



log(A fc ) = -— (- I ||V0|| 2 <M+ / k^dA + HA^+Cig), 

07T \Z JM JM / 



we have 
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Corollary 4.4 Let g and h be conformally equivalent and let 

(det(A h )\ 



M(h) = E(h) - 127rlog 



Then 1 9(h) = M{h) - M(g). 

So up to a constant our objective function is given by M(h). 



5 Questions 

Question 5.1 In |3|], Colin de Verdiere suggested that objective functions 
related to circle pattern problems might be related to the determinant of the 
Laplacian. One certainly would hope for a more direct relationship than that 
provided by corollary ^4.4| , raising the following natural question: Is there a 
disk pattern objective function, D^, with the property that E^(D/J limits to 
log(det(Aft)) when suitably normalized? 



Question 5.2 Corollary [4.4| also reveals some unexpected properties of I 9 . 
For example, there was no reason to expect that the roles of the g and h 
metric could be decoupled via I 9 {h) = M(h) — M(g); or even the mysterious 
implication that I 9 {h) = —I h (g). Furthermore corollary [L4| allows us to ex- 
tend this objective function consistently to the space of all negatively curved 
metrics. Notice that I 9 {h) can be computed and discretely interpreted when 
g and h are not in the same conformal class. With this observation we are 
left with the following natural question: Does corollary [4.4| hold among all 
metrics with negative curvature? Notice in particular this would provide a 
local formula for the log(det(A/ l )) outside a conformal class. 



Question 5.3 The idea of using disk patterns to explore uniformization can 



be traced back to Thurston Thurston's original idea was to approximate 
the Riemann mapping with a disk pattern solution. Thurston's idea was 
initially justified in [[13] and has been developed considerably since then, 
see for example [[7|. It would be interesting to implement the spirit of this 
approximation approach in this setting. For example it would be nice to 
answer the following question: If one takes a "dense" set of points on a 
Riemannian surface and forms the "empty disk pattern", can one measure 
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in a meaningful way how far the uniform surface produced by theorem 33 
is from the conformally equivalent uniform structure produced by theorem 
O? 



Question 5.4 This whole story carries over to the spherical case, with the 
exception of the convexity of the analogs of the objective functions H and I 9 . 
Can one use these same techniques to prove the corresponding uniformization 
results, even without convexity? 

Question 5.5 Can techniques like these be applied to find geometric struc- 
tures on 3-manifolds? 
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